Exotic leaky wave radiation from anisotropic epsilon near zero metamaterials 
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We investigate the emission of electromagnetic waves from biaxial subwavelength metamaterials. For tunable 
anisotropic structures that exhibit a vanishing dielectric response along a given axis, we find remarkable variation 
in the launch angles of energy associated with the emission of leaky wave radiation. We write closed form 
expressions for the energy transport velocity and corresponding radiation angle cp, defining the cone of radiation 
emission, both as a functions of frequency, and material and geometrical parameters. Full wave simulations 
exemplify the broad range of directivity that can be achieved in these structures. 
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Metamaterials are composite structures engineered with 
subwavelength components, and with the purpose of manipu- 
lating and directing electromagnetic (EM) radiation. Recently 
many practical applications have emerged, and structures fab- 
ricated related to cloaking, metamaterial perfect absorbers [1], 
and chirality [2, 3]. For metamaterials, the desired EM re- 
sponse to the incident electric (E) and magnetic (H) fields, 
typically involves tuning the permittivity, e, and permeability, 
/i in rather extraordinary ways. This includes double nega- 
tive index media (negative real parts of both e and /i), single 
negative index media (negative real part of e or ji), matched 
impedance zero index media [4, 5] (real part of e and ji is near 
zero), and epsilon near zero (ENZ) media (real part of e is 
near zero). Scenarios involving ENZ media in particular have 
gained prominence lately as useful components to radiative 
systems over a broad range of the EM spectrum [6-8]. 

In conjunction with ENZ developments, there have also 
been advances in infrared metamaterials, where thermal emit- 
ters [9], optical switches [10], and negative index metamate- 
rials [2, 11] have been fabricated. Due to the broad possibili- 
ties in sensing technologies, this EM band is of considerable 
importance. Smaller scale metamaterial devices can also of- 
fer more complex and interesting scenarios, including tunable 
devices [12], filters [13], and nanoantennas [14]. For larger 
scale ENZ metamaterials, high directivity of an incident beam 
has been demonstrated [6] . This can be scaled down and ex- 
tended to composites containing an array of nanowires, yield- 
ing a birefringent response with only one direction possessing 
ENZ properties [13]. A metamaterial grating can be designed 
to also have properties akin to ENZ media [15]. 

Often times, the structure being modeled is assumed 
isotropic. Although this offers simplifications, anisotropy is 
an inextricable feature of metamaterials that plays a crucial 
role in their EM response. For instance, at optical and in- 
frared frequencies, incorporating anisotropy into a thin planar 
(nonmagnetic) waveguide can result in behavior indicative of 
double negative index media [16]. Anisotropic metamaterial 
structures can now be created that contain elements that pos- 
sess extreme electric and magnetic responses to an incident 
beam. The inclusion of naturally anisotropic materials that 



are also frequency dispersive (e.g., liquid crystals), allows ad- 
ditional control in beam direction. It has also been shown that 
metamaterial structures requiring anisotropic permittivity and 
permeability can be created using tapered waveguides [17]. 
By assimilating anisotropic metamaterial leaky wave struc- 
tures within conventional radiative systems, the possibility ex- 
ists to further control the emission characteristics. 

Prompted by submicron experimental developments, and 
potential beam manipulation involving metamaterials with 
vanishing dielectric response, we investigate a planar 
anisotropic system with an ENZ response at near-ir frequen- 
cies along a given (longitudinal) direction. By "freezing" 
the phase in the longitudinal direction and tuning the elec- 
tric and magnetic responses in the transverse directions, we 
will demonstrate the ability to achieve remarkable emission 
control and directivity. When excited by a source, the direc- 
tion of energy flow can be due to the propagation of localized 
surface waves. There can also exist leaky waves, whereby the 
energy radiatively "leaks" from the structure while attenuating 
longitudinally. Indeed, there can be a complex interplay be- 
tween the different type of allowed modes whether radiated or 
guided, or some other mechanism involving material absorp- 
tion. Through a judicious choice of parameters, the admitted 
modes for the metamaterial can result in radiation launched 
within a narrow cone spanned by the outflow of energy flux. 

Some of the earliest works involving conventional leaky 
wave systems reported narrow beam width antennas with pre- 
scribed radiation angles [18], and forward/backward leaky 
wave propagation in planar multilayered structures [19]. In 
the microwave regime, photonic crystals [9, 20, 21] and trans- 
mission lines can also can serve as leaky wave antennas [22]. 
More recently, a leaky wave metamaterial antenna exhibited 
broad side scanning at a single frequency [22]. The leaky 
wave characteristics have also been studied for grounded sin- 
gle and double negative metamaterial slabs [23]. Directive 
emission in the microwave regime was demonstrated for mag- 
netic metamaterials in which one of the components of is 
small [24] . Nonmagnetic slabs can also yield varied beam di- 
rectivity [25]. 

To begin our investigation, a harmonic time dependence 
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exp(— icjt) for the TM fields is assumed. The planar structure 
contains a central biaxial anisotropic metamaterial of width 
2d sandwiched between the bulk superstrate and substrate, 
each of which can in general be anisotropic. The material 
in each region is assumed linear with a biaxial permittivity 
tensor, = ef^xx + c(^yy + cI^t^t^. Similarly, the bi- 
axial magnetic response is represented via iii = jj.f^x.-k + 
/i^^yy + /if^zz. The translational invariance in the y and z 
directions allows the magnetic field in the ith layer, H^, to 
be written = yh\{x)e^^^^~'^'^\ and the electric field as 
= [xe^^{x) + zel{x)\e'^^^-^^\ Here, 7 = /3 + m is 
the complex longitudinal propagation constant. We focus on 
wave propagation occurring in the positive x direction, and 
nonnegative P and a. Upon matching the tangential E and 
H fields at the boundary, we arrive at the general dispersion 
equation that governs the allowed modes for this structure, 

^r^±,2(erfe±,i+effe^,3)+ (1) 
[(ef )2^^,i^^,3 - ef ef ^1,2] tan(2d/c^,2) = 0, 

where the transverse wavevector in the superstrate (referred to 
as region 1), is, 

(2) 

For the metamaterial region (region 2), we write k^^2 = 
-^^k^ji^ — 7^62^/62^, and for the substrate (region 3), 
/c^,3 = ±^7^63^/63^ — k^ji^^el^. The choice of sign in re- 
gions 1 and 3 plays an important role in the determination of 
the physical nature of the type of mode solution that will arise. 
The two roots associated with k^^2, results in the same solu- 
tions to Eq. (1). The dispersion (Eq. (1)) is also obtained from 



the poles of the reflection coefficient for a plane wave incident 
from above on the structure. The transverse components of the 
E field in region 1 are, ef = -ik^^i/ (kQel^)Hie~^^^^^'^~'^\ 
ef = 7/(A:oef^)/if , and h\ = i^ie'-^^'^^^-^), where Hi is a 
constant coefficient. 

Next, to disentangle the evanescent and leaky wave fields, 
we separate the wavevector k±^i into its real and imagi- 
nary parts: k±^i = + ^^^), with g+ and q~ real. 
The k±^i, q~, and are in general related, depending on 
sgn(ef^a/3/ef^). For upward wave propagation (-\-x direc- 
tion), clearly we have q~^q~ > 0. It is also apparent that the 
parameter q~ represents the inverse length scale of wave in- 
crease along the transverse x-direction. We are mainly con- 
cerned with the k±^i that correspond to exponential wave 
increase in the transverse direction while decaying in 2;, a 
hallmark of leaky waves. Although leaky wave modes are 
not localized, they can be excited by a point or line source 
which gives rise to limited regions of space of EM wave 
amplitude increase before eventually decaying. When ex- 
plicitly decomposing k±^i into its real and imaginary parts, 
there is an intricate interdependence among 7, e^, and fii 
(for a ^ 0): q^ = 1/V2{V^^TW T Af'^. where 
A = efVef^(/32 - a^) - e^^ and B = 2a/3efVef^. 
We will see below that q^ is the root of interest in determining 
leaky wave emission for our structure. At this point the sur- 
rounding media can have frequency dispersion in e^, and fii, 
while the anisotropic metamaterial region can be dispersive 
and absorptive. 

We are ultimately interested in anisotropic metamaterials 
with an ENZ response along the axial direction (2;-axis). In the 
limit of vanishing e^^, and perfectly conducting ground plane, 
Eq. (1) can be solved analytically for the complex propagation 
constant, 7. The result is 
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The two possible roots correspond to distinct dispersion 
branches (seen below). There are, in all, four solutions, 7^, 
and —7^. The geometrical and material dependence con- 
tained in Eq. (3), determines the entire spectrum of the leaky 
radiation fields that may exist in our system. 

There are numerous quantities one can study in order to 
effectively characterize leaky wave emission. One physi- 
cally meaningful quantity is the energy transport velocity, 
vt, which is the velocity at which EM energy is transported 
through a medium [26, 27]. It is intuitively expressed as the 
ratio of the time-averaged Poynting vector, 6'avg, to the energy 
density, U: Vt = S^vg/U. Properly accounting for frequency 
dispersion that may be present, we can thus write the energy 



velocity for EM radiation emitted above the structure, 

c/(87r)Re[£:i x H^] 
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where the conventional definition [28] of U has been extended 
to include anisotropy. Inserting the calculated EM fields, we 
get the compact expression (assuming no dispersion in the su- 
perstrate). 



(5) 



The corresponding direction of energy outflow is straightfor- 
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FIG. 1: (Color online). The real and imaginary (a) parts of 
the complex propagation constant 7^, normalized by the vacuum 
wavevector ko at f = 280 THz (/i^^ = 1). The figures (a) and (c) 
are 3D global views depicting a and /3 as functions of and the 
thickness parameter d. Figures (b) and (d) represent the normalized 
f3 and a, respectively, as functions of ef"^ and for d = 0.01 |im (solid 
curve), d = 0.05 |im (dotted curve), and d = 0.1 |im (dashed curve). 

wardly extracted from the vector directionality in Eq. (5), 

^ = tan-n^), (6) 

which holds in the case of loss and frequency dispersion in 
the metamaterial. It is evident that Eq. (6) satisfies ^ {) 
as a 0, corresponding to the disappearance of the ra- 
diation cone and possible emergence of guided waves. In 
this limit, Vt = zuj/ which corresponds to the expected 
phase velocity, or velocity at which plane wavefronts travel 
along the z-direction. There is also angular symmetry, where 
(^{€2^) (p{—e2^), when fi2^ — ^2^- For high refrac- 
tive index media (ef^ or /i2^ 00), we moreover recover the 
expected result that (p tends toward broadside ((p = 0). 

Leaky wave emission from an anisotropic metamaterial in 
vacuum is characterized in Figs. 1 (a) and (c), where 3-D 
views of the normalized (3 (Re [7+]), and a (Im[7+]), are 
shown as functions of the transverse dielectric response, ef^, 
and thickness parameter, d (the width = 2d). In Fig. 1 (b) 
and (d), 2D slices depict the normalized (3 and a as functions 
of ef^. Only the positive root, 7+, is shown, corresponding 
to the leaky wave case of interest, with a > 0. The slight 
kinks in the curves (see Fig. 1(b)) are at points where the 
7~ solutions would emerge (for a < 0). Both panels on 
the left clearly demonstrate how (3 rises considerably with in- 
creasing lef^l. For subwavelength widths (kod <C 1), and 
to lowest order, the propagation constant varies linearly in 
ef^, as /3/ko ^ ef^/ {2kody/ef^el^). As the dielectric re- 
sponse ef^ vanishes, corresponding to an isotropic ENZ slab, 
we see from the figures (and Eq. (3)) that /3 ^ (long wave- 
length limit), and emission is subsequently perpendicular to 



FIG. 2: (Color online) Leaky wave launch angle, cp, as a func- 
tion of permittivity, (Figs, (a) and (b)) for eight differ- 
ent thicknesses in succession, starting with d — 1 jam (dot- 
ted orange curve), and subsequent values of d (in |im), equal- 
ing 1/5, 1/10, 1/20, 1/30, 1/40, 1/50, and 1/60. Other param- 
eters are as in Fig. 1. In (c) the emission angle is shown as a 
function of frequency for the same thicknesses in (a) and (b). In 
(d), the effects of geometrical variation are presented for ef"^ = 
0.001, 0.01, 0.05, 0.2, 0.4, 0.6, and 0.8. The curves with larger over- 
all (p correspond to smaller ef^ in succession. 

the interface (see below). It is also interesting that the im- 
portant parameter a characterizing leaky waves rapidly in- 
creases from zero at ef^ = and peaks at differing values, 
depending on the width of the emitting structure (Figs. 1 (c) 
and (d)), until eventually returning to zero at the two points, 
= 4[-2{kod)'^ ± ^(MF+^^MF]- This illustrates 
that a is spread over a greater range of ef^ for larger widths, 
but as previously discussed in conjunction with Fig. 1, a si- 
multaneously suffers a dramatic reduction. For small d/X, the 
extremum of Eq. (3), reveals that the strength of the a peaks, 
amax, are given by amax ^ 1/2 ± kod. 

Next, in Fig. 2, we show the angle, (p, which defines the 
radiation cone from the surface of the metamaterial structure, 
as functions of both ef^, frequency, and thickness parame- 
ter d. In panel (a), the variation in (p is shown over a broad 
range of ef^ for nonmagnetic media (/i^^ = 1), while panel 
(b) is for a metamaterial with vanishing /i^^, representative of 
a type of matched impedance [5]. The eight curves in Fig. 2 
(a) and (b) represent different widths, identified in the cap- 
tion. We see that for ef^ 0, we recover the isotropic result 
of nearly normal emission ((p ^ 90°), discussed and demon- 
strated in the millimeter regime [6]. This behavior can be un- 
derstood in our system, at least qualitatively, from a geomet- 
rical optics perspective and a generalization of Snell's Law 
for bianisotropic media [29]. When the magnetic response 
vanishes (Fig. 2 (b)), the emission angle becomes symmetric 
with respect to ef^, dropping from (p = 7r/2 for zero ef^, to 
broadside ((p = 0) when ef^ = ±4/co<i. Thus thinner widths 
result in more rapid beam variation as a function of ef^. In 
Fig. 2(c) we show how the emission angle varies as a func- 
tion of frequency, with the transverse response obeying the 
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FIG. 3: Normalized field profiles illustrating broad angular variation 
in beam emission. The arrows along the interfaces depict the Poynt- 
ing vector. The top left and right panels correspond to ^ = 0.05, 
and ef"^ = 0.05 + 0.02i, respectively. The the bottom left and right 
panels are for ef^ = 0.66 and ef^ = 0.66 + 0.02i, respectively. 
The metamaterial is subwavelength (d = 1/20 |im) and nonmag- 
netic = 1). Coordinates are given in units of ( x 10) |im. 



Drude form, ef ^ = 1 - ujI/{uj^ + iVuj). Here, ujp = (27r)120 
THz and T = 0, to isolate leaky wave effects. With in- 
creasing frequency, we observe similar trends found in the 
previous figures, where a larger dielectric response pulls the 
beam towards the metamaterial. In (d), a geometrical study 
illustrates how the emission angle varies with thickness: for 
^xx^yy ^ 1^ the emission angle rises abruptly with increased 
d, before leveling off at = tan~-^(Y^l/(e2^/i2^) — 1). 
Physically, as the slab increases in size, the complex prop- 
agation constant becomes purely real, 7 y^ef^/if^, and 
Y^T^^e|^/i|^- This is consistent with what was dis- 
cussed previously involving the depletion of a with d; for 
thick ENZ slabs, leaky wave radiation is replaced by conven- 
tional propagating modes. For fixed ef^, there is also a critical 
thickness, d*, below which no leaky waves are emitted, which 
by Eq. (3) is, d* = e|^/(4A:o vT^^ef>f ). These results are 
consistent with simulations from a commercial finite element 
software package [30]. In Fig. 3, we show the normalized \ H\ 
arising from a source excitation (at / = 280 THz) within the 
metamaterial for d = 1/20 |im. The left two panels are for 
F = 0, and the right two have absorption present. The full 
wave simulations agree with Fig. 2(a) (dashed green curve), 
where the leaky-wave energy outflow spans a broad angular 
range when < ef^ < 0.68. The right two panels exem- 
plify the robustness of this effect for moderate amounts of loss 
present in the metamaterial. 



In conclusion, we have demonstrated leaky wave radiation 
in subwavelength biaxial metamaterials with vanishing per- 
mittivity along the longitudinal direction. The leaky- wave ra- 
diation cone illustrated broad directionality through variations 
in the transverse EM response. By utilizing nanodeposition 
techniques, such structures can be fabricated by implementing 
an array of metallic nanowires embedded in a self-organized 
porous nanostructured material [13]. 

K.H. is supported in part by ONR and a grant of HFC re- 
sources as part of the DOD HPCMP. 



[1] X. Liu, et aL, Phys. Rev. Lett. , 104, 207403 (2010). 
[2] S. Zhang, etaL, Phys. Rev. Lett. , 102, 023901 (2009). 
[3] C. Menzel, et ai, Phys. Rev. Lett. , 104, 253902 (2010). 
[4] R.W. Ziolkowski, Phys. Rev. E 70, 046608 (2004). 
[5] V. C. Nguyen, et aL, Phys. Rev. Lett. 105, 233908 (2010). 
[6] S. Enoch, et al, Phys. Rev. Lett. 89, 213902 (2002). 
[7] M. Silveirinha and N. Engheta, Phys. Rev. Lett. 97, 157403 
(2006). 

[8] A. Alu, et al, Phys. Rev. B , 65, 155410 (2007). 

[9] M. Laroche, et al, Phys. Rev. Lett. 96, 123903 (2006). 
[10] N.-H. Shen, et al, Phys. Rev. Lett. , bf 106, 037403 (2011). 
[11] S. Zhang, etal, Phys. Rev. Lett. , 95, 137404 (2005). 
[12] N.-H. Shen, et al, Phys. Rev. Lett. 106, 037403 (2011). 
[13] L.V. Alekseyev, etal, Appl. Phys. Lett. 97, 131107 (2010). 
[14] X-X. Liu and A. Alu, Phys. Rev. B 82 144305 (2010). 
[15] V. Mocella, et al. Opt. Express 18, 25068 (2010). 
[16] V.A. Podolskiy and E.E. Narimanov, Phys. Rev. B, 71, 

201101(R) (2005). 
[17] 1. 1. Smolyaninov, etal, Phys. Rev. Lett. 102, 213901 (2009). 
[18] R. E. Collin and F. J. Zucker, Antenna Theory, Part 11, Mcgraw- 
Hill (1969). 

[19] T. Tamir and E Y. Kou, IEEE J. Quantum Electr. 22, 544 (1986). 
[20] E. Colak, et al. Optics Express 17, 9879 (2009). 
[21] A. Micco, et al, Phys. Rev. B 79, 075110 (2009). 
[22] S. Lim, IEEE trans, on Microw. theory and Tech. 52, 2678 
(2004). 

[23] P. Baccarelli, et al, IEEE Trans, on Microw. Th. and Tech., 53, 
32 (2005). 

[24] Y. Yuan, et al, Phys. Rev. A 77, 053821 (2008). 

[25] H. Liu and K. J. Webb, Phys. Rev. B 81, 201404(R) (2010). 

[26] R. Loudon, J. Phys A 3, 233 (1970). 

[27] R. Ruppin, Phys. Lett. A 299, 309 (2002). 

[28] L. D. Landau, E. M. Lifshitz, and L. P. Pitaevskii, Electro- 
dynamics of Continuous Media (Butterworth-Heinemann, Ox- 
ford, 1984), 2nd ed. 

[29] T. M. Grzegorczyk, et al, IEEE Trans. On Microw. Theory and 
Tech. 53, 1443 (2005). 

[30] COMSOL Multiphysics, http://www.comsol.com 



